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We extend the theory of Bose-Einstein condensation from Bogoliubov's weak-coupling regime to 
arbirarily strong coupHngs. 
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1. For (/)'^-theory in _D < 4 Euclidean dimensions with 
0(A'^)-synimetry, a powerful strong-coupling theory has 
been developed in 1998 It has been carried to 7th or- 
der in perturbation theory in D = 3 , and to 5th order 
in D = 4 — e dimensions The theory is an exten- 
sion of a variational approach to path integrals that was 
once set by R.P. Feynman and collaborator in 1989 
The extension to high orders is descibed in the textbook 
0. It is called Variational Perturbation Theory (VPT). 
Originally, the theory was designed to convert only the 
divergent perturbation expansions of quantum mechanics 
into exponentially fast convergent exressions @. In the 
papers [lj-l3| , it was extended from quantum mechanics to 
(/)^-theory with its anomalaus dimensions and produced 
all critical exponents. This is called quantum field theo- 
retic VPT. That theory is explained in the textbook 
and a recent review Isl. 



Surprisingly, this successful theory has not yet been 
applied to the presently so popular phenomena of Bose- 
Einstein condensation. These have so far mainly been fo- 
cused on the semiclassical treatments using the good- 
old Gross-Pitaevskii equations, or to the weak-coupling 
theory proposed many years ago by Bogoliubov (lo| . This 
is somewhat surprising since the subject is under intense 
study by many authors. So far, only the shift of the crit- 
ical temperature has been calculated to high orders [ll| . 
There are only a few exceptions. For instance, a simple 
extension of Bogoliubov's thepry to str ong couplings was 
proposed in [l^l and pursued further in [ij] . But that had 
an unpleasant feature that it needed two different chemi- 
cal potentials to maintain the long- wavelength properties 
of Nambu-Goldstone excitations required by the spon- 
taneously broken U(l)-symmetry in the condenste. For 
this reason it remained widely unnoticed. Another no- 
table exception is the theory in [l5[ which came closest 
to our approach, since it wa also based on a variational 
optimization of the energy. But by following Bogoliubov 
in identifying oq as ^/po from the outset, they ran into 
the notorious problem of violating the Nambu-Goldstone 
theorem. Another approach that comes close to ours is 
found in the paper [l6|. Here the main difference lies in 
the popular use the Hubbard-Stratonovic transformation 
(HST) to introduce a fluctuating collective pair field [17 1 . 
But, as pointed out in 24 1 and re-emphasized in this 
makes it impossible to calculate higher-order corrections 



The rules for applying VPT to nonrelativistic quantum 
field theories in 3+1 dimensions have been specified some 
time ago [isj . In this note we want to show how derive 
from them, to lowest order, the properties of the Bose- 
Einstein condensation at arbitrarily strong couplings. 

It must be mentioned that in the literature, there have 
been many attempts to treat the strong-coupling regime 
of various field theories for models with a large num- 
ber of identical field components (the so-called large- TV- 
models). This has first been done for the so-called spher- 
ical model [31, later the Gross- Neveu model [i^], and 
0(A^)-symmetric lyS^-models 2l|. In all these applica- 
tions, the leading large-A^ limit has been easily solved 
with the help of the HST trick of introducing a fluc- 
tuating fleld variable [13, [l^l for some dominant col- 
lective phenomenon (Collective Quantum Field Theory 
17|). This approach has, however the above-discussed 
problems of going to higher orders [25j , which are absent 
here. 

2. The Hamiltonian of the boson gas has a free term 



p 



(1) 



where Sp = p-^/2A/ are the single-particle energies and 
= Ep — /i the relevant energies in a grand-canonical en- 
semble. As usual, and Op are creation and annihilation 
operators defined by the canonical equal-time commuta- 



tors of the local fields ipix) 
interaction is 
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Instead of following Bogoliubov in treating the p = 
modes of the operators Op classically and identifying with 
the square-root of the condensate density pq, we intro- 
duce the field expectation = y^VSo/g as a varia- 



tional parameter , and rewrite ffint as Hi 
plus 
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(3) 



plus a fiuctuation Hamiltonian iJ;'^^. , which looks like ([2]) , 
except that the sum contains only nonzero-momentum 



2 



modes. Now we proceed according to the rules of VPT 
psj and introduce dummy variational parameter E and 
A via an auxiliary Hamiltonian 



trial ' 



in 



(4) 



leading a harmonic Hamiltonian 
V 
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^(ep - M + 2Eo)a|,ap 
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h.c. 



trial; 



(5) 



for which we have to calculate the energy prder by order 
in pcrturation theory considering 
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as the interaction Hamiltonian. The zeroth-order varia- 
tional energy is Wq = (^^o)' ^^'^ lowest-order correc- 
tion comes from the expectation value AiVF = (i/;™^' ). 
If the energy is calculate to all orders in H^^^ the result 
will be independent of the variational parameters Eqj 
and A, but the energy to any finite order will depend 
on it. The optimal values of the parameters are found 
by optimization (usually extremization), and the results 
converge exponentially fast as a function of the order 

A Bogoliubov transformation with as yet undeter- 
mined coefficients itp, Vp constrained by the condition 
1, produces a ground state with vacuum ex- 
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pectation values {a^Up/ 
that 



Vp and 



Wo 



V . 



9 2g " 
J2{[£p-^^-^ + 2So] fp + (Eo - A)upVp}. 



(7) 
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The first-order variational energy Wi contains, in ad- 
dition, the expectation value (i/;™^''). Of this, the first 
part, A(i o)W^ = {H['^^), is found immediately with the 
help of the standard commutation rules as a sum of three 
pair terms 



p+q"p'- 



p'flp; 



+ («p+qap)(aL Op,) + (aL ap-)(a!,, Op). (8) 



so that 



(9) 



The second part 
value 

A(i,i)W^ 



p.p'#o 

-iJtriai) adds to tlus the expectation 



PT^O 



Swp + AupWp) 



(10) 



In order to fix the average total number of particles 
iV, we differentiate Wi=Wq + /^(ifi]W + A(i_i)M^ with 
respect to — /i and set the result equal to N to find the 
density p = N /V as 



P-y + E-P- 

^ p#0 



(11) 



The momentum sum is the density of particles outside 
the condensate, the uncondensed density 



Pu = EK«p) = -^E^p 

Pt^O p 

implying that Eq /g is the condensate density po : 
Eo 
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= Po = P - PVL- 



(12) 



(13) 



Now we extremize Wi with respect to the variational 
parameter Eq which yields the equation 

^-E(2Wp + UpWp) = 2pu+E UpVp = 2p^ + 5. (14) 
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We are now able to fix the size of the Bogoliubov co- 
efficients Up and Vp. The original way of doing this is 
algebraic, based on the elimination of the off-diagonal el- 
ements of the transformed Hamiltonian operator. In the 
framework of our variational approach it is more natural 
to use the equivalent procedure of extremizing the energy 
expectation Wq with respect to Up and Vp under the con- 
straint Up — lip = 1, so that dup/dvp = Vp/up. Varying 
Wo, we obtain for each nonzero momentum the equation 

2(ep-Ai+2Eo-E)up+(Eo-A)(up + u2/up) =0. (15) 
In order to solve this we introduce the constant 

E=-Ai+2Eo-E =-/i + 2.g(p- pu) - E, (16) 



the right-hand side emerging after using (|T2|) and (|13p . 
We further introduce the constant 

A = En - A = rS. (17) 

Then we rewrite (jisp in the simple form 

2(ep + S)t;p + rE(up + i;2/up) =0, (18) 

which is solved for all p by the Bogoliubov transformation 
coefficients 



Ep + E 



Sr. 



, (19) 
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with UpVp = A/2£p, and the quasiparticle energies 



r2E2 



(20) 
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Having determined the Bogoliubov coefficients, we can 
calculate the above momentum sums in Eqs. (|12[) and 
dH]). We begin with the uncondensed particle density 
([T^ . Inserting ([T^ . it becomes 



_ 1 2 _ 1 /■ 



- 1 



(21) 



The integral is easily done if wc set |p| = hkj^ k with 

(22) 



fcg = V2M^/h, so that wc find 
where 



^2 + 1 



v/(k2 + 1)2 



-1^=^/W. (23) 



with /pi'' — 1. The second momentum sum in Eq. (|14p 
reads, after inserting p^ . 



(Pp 1 

2 J {2TThY&p' 

(24) 

In contrast to pT|) . this is a divergent quantity. As a 
consequence of the renormalizability of the theory, the 
divergence can be removed by absorbing it into the in- 
verse coupling constant of the model defined by 

1 _ 1 1^1 1 [ d^p I 



p#0 



(2^)3 2ep' 



The renormalized coupling is finite and measurable in 
two-body scattering as an s-wave scattering length: gn — 
2TTh'^as/M. Thus we introduce the finite renormalized 
quantity 



1 rH f d'^ 

fe = 7L«p-p--Ti (2^ 



rH f d^p f I 1 



(27rft)3 V'^^P £ 



(26) 



and write 5 = 5r + f^div, where the divergence is the 
momentum sum 



rS f d^p 1 



V ^2e 



2 J {2T:hf Ep 



If we denote this by — S/Fu, we have 

rE 

5 = 5r + (Jdiv = 5r - — . 

Vv 

Inserting this together with ([T^ into (|14p . we find 

^ - So 
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= 2pu + Sr + S, 



div ■ 



(27) 



(28) 



(29) 



Recalling p^ . this implies 



Po + 2pu + Sr + 6dw ^ p + Pm + 6r + (5div (30) 



If we evaluate the momentum sum (|26p in the same 
way as ((2T|) . it yields 



where is given by the integral 

1 



(31) 



-r I dK K 
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y =x/2/W,(32) 



with /]^^ = 1. 

We continue the discussion with Eq. (jl6p . which we 
rewrite using (|30p as 



9 9 Vv 



(33) 



As before in Eqs. ([25]), and ^6^, the last, divergent 
term can be absorbed into the first by renormalizing the 
coupling constant, so that wc obtain 



= P - oPu - Or . 

9R 9 



(34) 



Finally, we calculate the total variational energy Wi. 
Inserting the Bogoliubov coefficients ([T^ into Wq of 
Eq. ([7]) and adding the action energies A(i o)Ty+A(i W 
of and © and PH)) . wc have 



Wi 



— /iSo + t^Sq + ?«o(S, r) 



.g- 2g 



4Fu 



where wo(E) is the convergent momentum sum 



«;o(E,r) 



5I: 

p#0 



2e. 



This is evaluated as in (PT|) to 



(35) 



(36) 



(37) 



where 





8\/2 (r) 

15 ■'^ ' 



^(^2 + l)2_^2_^2_^^J_ 



(38) 



with f^' = I. If we rename all Ie/^t^"^ to J^;, the energy 



(39) 



becomes 

V V ^ V - ^ -(r) 

-H^fc|(2/M2+A/l^»2)^^^3/2^^jM+Ajf ). 

The expression is renormalized most simply using dimen- 
sional regularization, that allows us to use Veltman's rule 
to set 1/v = 0. 
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We now prepared to extremize the variational energy 
Eq. ((39)) with respect to E and A. We insert Yi/g = 
P ^ ^ E/p from ([M)) and IS./ g = p — — K/ g 
from pT]) and vary Vl^i in (5S and (5 A. This yields the 
equations 

(S/.g - p(-))5n + (A/5 - p(+))^i2 = 0, (40) 
{t/g - p(-))52i + (A/.g - p(+))522 = 0, (41) 

where pu'^' = P — Pu ± (^fl and 

Sii = i^;-]dkyd% s^2^{ky^i)di'f]idT, (42) 
= i^;Jdki/d% S22 = (fc|/E)a/i'-V9r. (43) 

These equations are solved for r = 1 with 



— = — = P - Pu - Or. 

9 9 



(44) 



The reason is simply that both IpJ and Jj behave near 
r = 1 like (1 — r)"^/^ so that derivative ar r = 1 vanishes 
and both equations (j43|) give (|44|) . 

The solution of these equations is ?' = 1, thus guaran- 
teeing the Nambu-Goldstone nature of the quansiparticle 
energies ([20]) . 

3. To extract experimental consequences it is useful 
to re-express all equations in a dimensionlcss form by 
introducing the reduced variables 



S 

s = — , 



(45) 



where = f? j^Ma} is the natural energy scale of the 
system. We also introduce the reduced s-wave scattering 
length 



as = OTT- 



(46) 



in terms of which the renormalized coupling constant is 

471^2 



9R 



M 



-as = S-KEaa as = Saa Os, (47) 



while 



a gu G'Ka'^as a^^Us 
and the second-sound velocity reads 



Is Pa h 

C = \l — Vn, V„ = = . 

2 ' M aM 



(48) 



(49) 



Let us also define the reduced quantity 6 = s^^'^Is and 
Pu = s3/%^. (50) 



In terms of these we calulate the reduced variational en- 
ergy wi = Wi/NSa from Eq. (|39l) for r = 1: 



/- \2 ^ 



Wi = -as(l + pu + ^)(l-Pu) + y(l-/5u)"' + ^-^£; 



y (2p2 +5'^)+ Gsia^Pu + (taS), 



(51) 



where crs = 1 — Sp^ — S — s/as from (|33p and cta = 
l—p^—s/cis from (|17p . Inserting here pu and 5, and going 
from the grand-canonical to the true proper energies by 
adding fiN to W\ forming — W\ -I- pV p^ we obtain 
the reduced energy 



4\/2 
157r2' 



,3/2 



- 3/2 , 1-4 



The relation between s and dg is from (|44p 



leading to expansion 



-5/2 



157r 



2 s 



727r' 



4"'s 



(52) 



(53) 



(54) 



Note that in the strong-coupling limit, s ^ s^'^ ~ 
(37r^/-\/2)^^'^, the maximal depletion is pu = 1/4, and 
the energy behaves like w\ — _B -|~ Ad^, with B = 
-4 X 21/12/(5 X 3i/*^7ri/3) w -0.48, and A = 1/2 + 
l/24^/27r2 + 2^/'^/\/iiT s a 0.72 . The sound velocity at 
infinite coupling is c = ^ s^'^ /2va ■ 

Let us now study the temperature dependence of our 
results. For this we introduce the temperature-dependent 
version of the integral at r = 1, where we omit the 
trivial superscript, to find 



Pu(i) =fc|^Pu(^)/47^^ fcs = V2MY./h = ^/^/a. (55) 
where Ip^ (t) is defined by the integral 



dKK 



ctin) - 1 



(56) 



where Ct(K) = coth \^\/{k'^ + 1)^ — l/2tj , and t is the re- 
duced temperature 



(57) 



t = keT/e^, = V2ME. 



To find the s at any temperature we need Eq. (j3ip for 
T 7^ 0, where it reads 



with 



5R{t) ^klh{t)/AT:\ 



(58) 



(59) 
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to find 



^ 47r2 ■ 



(60) 



For the temperature dependence of tlie energy in 
Eqs. (p6)) . and (|37| . we calculate 



Jo 



(61) 



The phase transition lies at the temperature where 
Pu = P- Along the transition, we find from Eq. (j53p 
the relation between s and a. 



3 /as = 



(62) 



For weak couplings, all expressions can be calculated 
analytically. The calculation is somewhat subtle since 
the small s -region of the integral cannot simply be ob- 
tained by expanding the integrand in powers of s, the 
first correction going like y/s. To see this, we must 
proceed as in the derivation of the Robinson expan- 
sion of the Bose-Einstein integral function [2^, writing 
Pu/p = s^/^l2/4:n^ + Apu/p, with the second term being 
the integral 



1 

4^ 



,(63) 



where f° = T^/Ta = 47rC(3/2)-2/3 ig the reduced critical 
temperature of the free Bose gas, and r is the ratio T/T^. 
The integral can be done immediately for s = and yields 
the well-known result 



p s^O p 

For small s, there is an additional subtracted term 



(64) 



Apu' 



1 

4^ 



dn 
2 



.2 



- 

2 



1 



(65) 



The first term takes its leading small-s behavior from the 
linear Nambu-Goldstone momentum behavior of second 
sound, becoming 



Apu' 2Tfo 



47r2 



1 



2sk2 



which is equal to — (tT" /47r) ■\/s/2. Thus we obtain for 
small s the leading terms of the uncondensed particle 
density 



P^i ^ 3/2 
p 3 • 47r2 



3/2 



' ^ c 

An 



(67) 



The last term is dominant for small s. Its negative sign 
has an interesting effect upon the phase diagram observed 
in earlier publications, that for small coupling constant, 
the critical temperature increases above the free Bose gas 



value to Tc = TcT° with 



--^+34^\/2+^(^)- 



(68) 



A similar limit square-root limit appears in the relation 
(f58| for Sji which becomes 



^ _ 3/2 ^sjt) _ 3/2 



p 47r2 47r2 47r2 

where the last term is equal to 



= s^,-y^ + ^s^/^Ahi{t), (69) 



47r2 



dn K 



1 



(70) 



To lowest order in r this yields 

so that we find 

^_^3/2/iW _ ^3/2^2 _rT;« /7 
p ^ 4^ ~ 4^2 4^ 

According to (|62p . this is equal to —s/ag, implying for 
small s, where the s^^^-tcrm is dominant, the relation 
between s and as along the phase transition line 



Inserting this into ([68)) . we obtain 



3 (47r): 



T&s + 



which becomes with f ° = [C(3/2)]-2/347r: 



Tc (Is 

TO " + 3C(3/2)4/3 



■ • = 1 



(73) 



(74) 



(75) 



where the constant is C w 87r/37rC(3/2)'^/3 w 2.33. This 
shows the initial increase of the critical temperature for 
small repulsion between the bosons discussed in [27 1. 
Numerically, the prcfactor C of the linear term is twice 

(66) 

as big as the value Cyp^^ — 0.93±0.13 predicted from 5 



loop variational perturbation theory in Ref. ll[ , and 83% 
larger than the value C = 1.27 ± 0.11 from its extension 
to seven-loops [l2j. It is the same as the value obtained 
from a large- calculation [2^ . 

For strong coupling, there is no subtlety and the phase 
transition can be extracted from a numerical plot of 
the location where ([50]) is equal to unity. Some plots 
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Figure 1: Transition temperature r = T/T^ as a function of 
a) s = Ti/sa, and b) of of the reduced s-wave scattering length 
as. 




Figure 2: a) Uncondensed density pu as a function of a) the 
quasigap parameter s and b) the reduced s-wave scattering 
length 0,1, at various reduced temperatures T /T^ = t — Q, 
0.25, 0.5. Note that the curves continue smoothly to negative 
l/cLs- 



are shown in Fig. [TJ The limit of infinitely strong cou- 
pling is found from the vanishing of the right-hand side 
of Eq. ([55]) . In particular we can easily calculate the 
temperature dependence of the second-sound velocity 

sftva as a function of temperature using Eq. (|53p . 

The behavior of Eqs. ([53]) and is shown in Figs. [21 
and [3l An interesting feature of Fig. [IJj) and Figs. [3] 
is that the curves continue smoothly beyond the strong- 
coupling limit l/cis = to negative values. This should 
be observable in experiments through a Feshbach reso- 
nance. 

4. All properties of the strongly interacting Bose gas 
determined by the above theory can be calculated in the 
presence of superflow of velocity f, by simply generalizing 
the integrals (f56| and (f59|) for /p„(t) and /^(i) to /p„(t, v) 
and /^(t, these being defined by interchanging in each 
integrand the terms (K^-t-1) by (k^ + 1 + 2ki//s^/^), where 
V is the reduced velocity of the gas v = v/va- From the 
associated second derivative of the energy we can easily 
find the superfluid density ps as a function of the velocity. 

There is no problem to drive the accuracy to any de- 
sired level, with exponentially fast convergence, as was 
demonstrated by the calulation of critical exponents in 
all Euclidean cp'^ theories with N components in D di- 
mensions [3l . The procedural rules were explained in the 
paper [25|. We merely have to calculate higher-order di- 
agrams using the harmonic Hamiltonian ([5]) as the free 
theory that determines the Feynman diagrams, and © 
as the interaction Hamiltonian that determines the ver- 
tices. At any given order, the results are optimized in 
the variational parameters T,q,T,, and A. The theory is 
renormalizable, so that all divergencies can be absorbed 
in a redefintion of the parameters of the orginal action, 
order by order. This is the essential advantage of the 



Pu/P 




Figure 3: a) Relation between the reduced s-wave scattering 
length cLs and quasigap parameter s = E/sa determining the 
second-sound velocity via c = ^JsJ^Va, at various reduced 
temperatures r = T jT'^ = 0, 0.5, 1, 1.5, and 2. b) Den- 
sity of uncondensed fluid pu as a function of s/fis at reduced 
temperatures r = 0, 0.25, 0.5, 0.75, and 1. 

present theory over any previous strong-coupling scheme 
published so far in the literature, in particular over those 
based on Hubbard-Stratonovic transformations of the in- 
teraction, which are applicable only in some large- A'^ limit 
as explained in [25j, and for which no higher- loop calcu- 
lations are renormahzable. 

Our results can be made much more reliable in the 
E ^ -regime by calculating the contribution of the still- 
missing second two- loop diagram.^ 
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